Uniform semiclassical approximations for umbilic bifurcation catastrophes 
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Gutzwiller's trace formula for the semiclassical density of 
states diverges at the bifurcation points of periodic orbits and 
has to be replaced with uniform semiclassical approximations. 
We present a method to derive these expressions from the 
standard representations of the elementary catastrophes and 
to directly relate the uniform solutions to classical periodic 
orbit parameters, thereby circumventing the numerical appli- 
cation of normal form theory. The technique allows an easy 
handling of ungeneric bifurcations with corank 2 such as the 
umbilic catastrophes and is demonstrated on a hyperbolic um- 
bilic in the diamagnetic Kepler problem. 

PACS numbers: 05.45.+b, 03.65.Sq, 32.60.+i 



Gutzwiller's periodic orbit theory ||,|^] has become the 
key for the semiclassical interpretation of quantum sys- 
tems with underlying chaotic classical dynamics. The 
contributions of isolated periodic orbits to the periodic 
orbit sum are given as 
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with T po , Sp , M po , and /x po the orbital period, classical 
action, stability matrix, and Maslov index, respectively. 
However, Eq. ^ diverges at the bifurcation points of pe- 
riodic orbits where orbits are not isolated, and (^) must 
be replaced with uniform approximations given in terms 
of diffraction catastrophe integrals ^jji). The study of 
bifurcations and uniform approximations is of fundamen- 
tal importance to the complete understanding and semi- 
classical treatment of systems with mixed regular-chaotic 
classical dynamics |5| . The derivation of the uniform so- 
lutions based on a canonical transformation of the coor- 
dinates and momenta to normal form coordinates [^|J^] 
and the construction of the diffraction integrals in terms 
of the new variables is usually a lengthy and nontriv- 
ial task P-f|[|, especially in the neighborhood of bifur- 
cations of codimension K > 2, and for catastrophes of 
corank 2 such as the umbilics. A simple scheme would 
be desirable to construct uniform approximations from 
classical periodic orbits and to relate the parameters of 
catastrophe diffraction integrals directly to the periodic 
orbit parameters, such as the classical action S and the 
eigenvalues of the stability matrix M. 

In this Letter we want to demonstrate that in prac- 
tical applications the derivation of uniform semiclassi- 
cal approximations can be considerably simplified espe- 
cially for ungeneric bifurcations of codimension K > 2 
and catastrophes of corank 2 when starting directly from 



the standard representation of the elementary catastro- 
phes • We illustrate our method by way of example 
of the diamagnetic Kepler problem Jll|,|l2| given by the 
Hamiltonian [in atomic units, 7 = 5/(2.35 x 10 5 T), and 
L z = m% = 0], 



h = -p 2 - - + - 7 y , 

2 y r 8 ' H ' 



(2) 



which is a scaling system, with w = 1/K e s = 7~ 1 / 3 the 
scaling parameter and E = i?7~ 2 / 3 the scaled energy. 
The classical dynamics is determined by the scaled en- 
ergy E but does not depend on w. We investigate the 
birth of four periodic orbits through two close-by bifurca- 
tions near the scaled energy E sa —0.096 where we search 
for both real and complex "ghost" orbits f|,[l3|. For the 
nomenclature of the real orbits we adopt the symbolic 
code of Ref. @. At scaled energy E { b 1} = -0.09689, the 
two orbits 00+- and +++ are born in a tangent bifur- 
cation. At energies E < E^\ a prebifurcation ghost 
orbit and its complex conjugate exist in the complex 
continuation of the phase space. Orbit 00+- is born 

unstable, and turns stable at the slightly higher energy 
~ (2) 

El = —0.09451. This is the bifurcation point of two ad- 
ditional orbits, 0-+ — and its time reversal +, which 

also have ghost orbits as predecessors. The graphs of the 
real orbits at energy E = are shown as insets in Fig. [I], 
and the classical periodic orbit parameters are presented 
as solid lines in Figs, [l] and ||. Fig. [I] shows the difference 
in scaled action, AS = AS/(2ttw), between the orbits. 

The action of orbit 0-+ — (or its time reversal +), 

which is real also for its prebifurcation ghost orbits, has 
been taken as the reference action. 

At this point the usual procedure would be to investi- 
gate the classical dynamics in the neighborhood of the pe- 
riodic orbits by numerical application of normal form the- 
ory HJtJ . The representation of the dynamics in normal 
form coordinates would finally lead to the correct type 
of the catastrophe diffraction integral related to the uni- 
form semiclassical approximation with numerically well 
determined coefficients. However, the numerical proce- 
dure of local canonical transformations to normal form 
coordinates, e.g., by means of local Fourier- Taylor series 
expansions with numerically obtained coefficients is 
rather lengthy and tedious especially for bifurcations re- 
lated to catastrophes of higher codimension or corank. 
The main result of this Letter is to demonstrate that 
there is a shortcut to the usual procedure which allows 
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FIG. 1. Difference AS* between the classical action of 
the four periodic orbits involved in the bifurcations. Dashed 
lines: Analytical fits related to the hyperbolic umbilic catas- 
trophe. Inset: Graphs of periodic orbits 0-+ — and its time 

reversal + (solid line), 00+- (dashed line), and +++ 

(dashed-dotted line) drawn in semiparabolical coordinates 
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circumventing the numerical application of normal form 
theory. By our new method, an easy construction of 
uniform scmiclassical approximations for ungeneric types 
of catastrophes, e.g. the umbilics, becomes feasible for the 
first time. 

Choosing the elementary catastrophe diffraction inte- 
grals as the ansatz for the uniform semiclassical approxi- 
mation, we must be able to identify the stationary points 
of the exponents with the periodic orbits, i.e., in our ex- 
ample four stationary points must exist. From the seven 
"elementary catastrophes" of Refs. ||,|lCj this is the case 
only for the swallowtail and the elliptic and hyperbolic 
umbilic. The correct choice in our example turns out to 
be the hyperbolic umbilic catastrophe, which is of im- 
portance, e.g., for uniform S matrix approximations in 
semiclassical scattering theory |l5| . It is a corank 2 catas- 
trophe, i.e., the diffraction integral is two-dimensional, 



dp 



dqe 
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with 



®(p, q; x, y) =p 3 + q 3 + y{p + qf + x[p + q) 



(3) 



(4) 



FIG. 2. Same as Fig. 1 but for the determinant det(M — 7) 
of the periodic orbits. 

For our convenience the function ^(p, q; x, y) slightly dif- 
fers from the standard polynomial of the hyperbolic um- 
bilic given in Ref. but the diffraction integral (|^) 
can be easily transformed to the standard representation. 
The four stationary points of the integral (^|) are readily 
obtained from the condition V<f> = as 
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= ±y/-x/3 => $(p , 9o; x,y) = 



(5) 



and 



2 + / 4 2 X 
PQ = qo = --y±^-y2-- 



.3/2 



t>( Pa ,q ;x,y) = ^y Qy 2 - xj t4 Qy 2 - |j . (6) 

The function Q(po, qo] x, y) must now be adapted to the 
classical action of the four periodic orbits, i.e., AS = 
2irwAS « &(po, qo', x, y), which is well fulfilled for 



x = aw 2/3 (E -E, 



=,(2) 
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and constants a = —5.415, b = 0.09665, as can be seen 
from the dashed lines in Fig. [l| Note that the agreement 
holds for both the real and complex ghost orbits. 

The next step to obtain the uniform solution is to cal- 
culate the diffraction integral (H) within the stationary 
phase approximation. For E > EYJ there are four real 
stationary points (po, qo) (see Eqs. [s] and |^) , and after ex- 
panding ^(p, q; x, y) around the stationary points up to 
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second order in p and q, the diffraction integral becomes 
the sum of Fresnel integrals, viz. 



2n 



E 



(4y 2 - 3x) T 2yx/^y 2 ~ 3a; 

(8) 



The terms of Eq. || can now be compared to the standard 
periodic orbit contributions (|]) of Gutzwiller's trace for- 
mula. In our example the first term is related to the 
orbit 0-+ — (with a multiplicity factor of 2 for its time 

reversal +), and the other two terms are related to 

the orbits 00+- and +++ for the upper and lower sign, 

respectively. The phase shift in the numerators describe 
the differences of the action AS and of the Maslov index 
Afj, = =pl relative to the reference orbit 0-+ — . The de- 
nominators are, up to a factor cw 1 / 3 , with c = 0.1034, 
the square root of | det(M — J)|, with M the stability ma- 
trix. Fig. ^| presents the comparison for the determinants 
obtained from classical periodic orbit calculations (solid 
lines) and from Eqs. and |§| (dashed lines). The agree- 
ment is very good for both the real and complex ghost 
orbits, similar to the agreement found for AS in Fig. [l]. 
The constant c introduced above determines the normal- 
ization of the uniform semiclassical approximation for the 
hyperbolic umbilic bifurcation which is finally obtained 
as 
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with To, So, and /xo denoting the orbital period, action 
and Maslov index of the reference orbit 0-+ — , and con- 
stants a, b, and c as given above. Note that all param- 
eters are readily determined by classical periodic orbit 
calculations. 

The comparison between the conventional semiclassi- 
cal trace formula (Q) for isolated returning orbits and 
the uniform approximation (g) for the hyperbolic um- 
bilic catastrophe is presented in Fig. || at the magnetic 
field strengths 7 = 10" 7 , 7 = 10" 8 , and 7 = 10" 9 , . 
For graphical purposes we suppress the highly oscillatory 
part resulting from the function exp[i(Sb/7i — f Mo)] by 
plotting the absolute value of A(E, w) instead of the real 
part. The dashed line in Fig. |] is the superposition of 
the isolated periodic orbit contributions from the four or- 
bits involved in the bifurcations. The modulations of the 

amplitude are caused by the constructive and destructive 

~ (2) 

interference of the real orbits at energies E > El ' and 
are most pronounced at low magnetic filed strength (see 
Fig. |]c). The amplitude diverges at the two bifurcation 
points. For the calculation of the uniform approximation 
(^) we numerically evaluated the catastrophe diffraction 
integral (0) using a more simple and direct technique as 
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FIG. 3. Semiclassical amplitudes (absolute values) for 
magnetic field strength (a) 7 = 10~ r , (b) 7 = 10~ 8 , and 
(c) 7 = 10 -9 in units of the time period To. Dashed 
line: Amplitudes of the standard semiclassical trace formula. 
Dashed-dotted line: Ghost orbit contribution. Solid line: Uni- 
form approximation of the hyperbolic umbilic catastrophe. 

described in ]jq| . Details of our method which is based 
on Taylor series expansions will be given elsewhere [n7|. 
The solid line in Fig. is the uniform approximation ^). 
It does not diverge at the bifurcation points but decreases 
exponentially at energies E < E^\ At these energies no 
real orbits exist and the amplitude in the standard formu- 
lation would be zero when only real orbits are considered. 
However, the exponential tail of the uniform approxima- 
tion (H) is well reproduced by a ghost orbit HQ with 
positive imaginary part of the complex action. As can 
be shown, the asymptotic expansion of the diffraction in- 
tegral (||) has, for x 3> 0, exactly the form of Eq. |] but 
with complex action S and determinant det(M — I) Jl7| . 
The ghost orbit contribution is shown as dashed-dotted 
line in Fig. ||. 

To verify the hyperbolic umbilic catastrophe in quan- 
tum spectra we diagonalized the Hamiltonian (||) in a 
complete basis set (for computational details see, e.g., 
JbSl) at constant scaled energy E = —0.1, which is 
slightly below the bifurcation energies, and calculated 
9715 eigenvalues w n for the scaling parameter in the re- 
gion w < 140. The scaled spectrum was analyzed by the 
high resolution method of Ref. i.e., the density of 
states g(w) = ^ n S(w — w n ) was fitted by application of 
the harmonic inversion technique to the functional form 
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of the semiclassical trace formula 

e {w)=Y d Me- 2 ™~ s * w (10) 

k 

with complex parameters and Sk- For isolated return- 
ing orbits these parameters, obtained from the quantum 
spectra, can directly be compared to the periodic orbit 
parameters of the classical calculations JlTJ] . The part of 
the complex action plane which is of interest for the hy- 
perbolic umbilic catastrophe discussed above is presented 
in Fig. ||. The two solid peaks mark the positions Sk and 
the absolute values of amplitudes \Ak \ obtained from the 
quantum spectrum. However, there is only one classical 
ghost orbit which is of physical relevance (dashed-dotted 
peak in Fig. ^). The position of that peak is in good 
agreement with the quantum result but the amplitude is 
enhanced, as is expected for isolated periodic orbit con- 
tributions near bifurcations (see Fig. y). For comparison 
we have analyzed the uniform approximation (^J) at con- 
stant scaled energy E = —0.1 and in the same range 
< w < 140 by applying the harmonic inversion tech- 
nique of Ref. [jig}] . The results for the uniform approx- 
imation are presented as dashed peaks in Fig. |J. The 
two peaks agree well with the quantum results for both 
the complex actions and amplitudes. The enhancement 
of the ghost orbit peak and the additional non-classical 
peak observed in the quantum spectrum are therefore 
clearly identified as artifacts of the bifurcation, i.e., the 
hyperbolic umbilic catastrophe. 
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FIG. 4. High resolution recurrence spectra at scaled en- 
ergy E — —0.1. Solid peaks: Part of the quantum recurrence 
spectra. Dashed-dotted peak: Classical ghost orbit contribu- 
tion. Dashed peaks: Uniform approximation of the hyperbolic 
umbilic catastrophe. 

In conclusion, we have presented a simple method to 
construct uniform approximations for the semiclassical 
density of states and to relate the parameters of catas- 
trophe diffraction integrals directly to periodic orbit pa- 
rameters such as classical action and eigenvalues of the 
stability matrix at energies near the bifurcation. The 
method is a shortcut to the conventional procedure, i.e., 
it circumvents the analysis of classical dynamics in the 
neighborhood of periodic orbits by numerical application 



of normal form theory, and therefore allows, for the first 
time, an easy handling of ungeneric bifurcations of sev- 
eral orbits related to catastrophes of higher codimension 
and corank. This has been demonstrated by way of ex- 
ample of a hyperbolic umbilic catastrophe in the diamag- 
netic Kepler problem, but evidently the method may be 
applied to other systems and catastrophe types as well. 
The technique will be useful for the semiclassical quan- 
tization of systems with mixed regular-chaotic classical 
dynamics, e.g., in combination with the method of har- 
monic inversion which has been successfully applied to 
systems with complete hyperbolic dynamics |2(J . 
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